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Abstract 

Weakly nonlinear plane waves are considered in hyperelastic crystals. Evolu- 
tion equations are derived at a quadratically nonlinear level for the amplitudes of 
quasi-longitudinal and quasi-transverse waves propagating in arbitrary anisotropic 
media. The form of the equations obtained depends upon the direction of propa- 
gation relative to the crystal axes. A single equation is found for all propagation 
directions for quasi-longitudinal waves, but a pair of coupled equations occurs for 
quasi-transverse waves propagating along directions of degeneracy, or acoustic axes. 
The coupled equations involve four material parameters but they simplify if the 
wave propagates along an axis of material symmetry. Thus, only two parameters 
arise for propagation along an axis of two-fold symmetry, and one for a three-fold 
axis. The transverse wave equations decouple if the axis is four-fold or higher. In 
the absence of a symmetry axis it is possible that the evolution equations of the 
quasi-transverse waves decouple if the third order elastic moduli satisfy a certain 
identity. The theoretical results are illustrated with explicit examples. 



1 Introduction 



We characterize and analyze degenerate weakly nonlinear elastic waves in crystals. The 
degeneracy considered here arises from the existence of acoustic axes in elastic materials. 
Acoustic axes are directions for which the phase velocities of at least two waves coincide. 
For classical elasticity this phenomenon typically occurs for transverse or quasi-transverse 
waves. In the mathematical literature the coincidence of wave speeds is called a loss of 
strict hyperbolicit y. Conditions for the existence of acoustic a xes were f irst d erived by 
KhatkevichI (1l962l ). and a useful re view of the topic is given by iFedorov f lQGSh. Recen t 
development s can be found in e.g. iBoulanger and HayesI f 1998l )r iMozhaev et all f 200ll ): 
NorrisI ( 20041 ). Analysis an d properties of nonlinear e lastic waves propagating alo ng acous- 
tic axes were discussed in IShuvalov and Radowicd f 200ll ). see also the book of iLyamov 
( 119831 ). 

The existence of acoustic axes, i.e. the loss of strict hyperbolicity, is typically accom- 
panied by the local loss of genuine nonlinearity, that is, vanishing of the scalar product 
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of the gradient of the phase velocity with the corresponding polarization vector, all eval- 
uated at the origin. The local loss of genuine nonlinearity implies the presence of weaker 
than quadratic (e.g. cubic) nonlinearities in the decoupled evolution equations for de- 
generate weakly nonlinear (quasi-)transverse waves. This happens e.g. for shear elastic 
plane waves in an isotropic mater ial as well as in a cubic crystal e.g. for [1 0] or [1 1 0] 
directions (see iDomanskil (j2000al )). However, coupled quadratically nonlinear evolution 
equations do occur for tran sverse elastic waves. Although this is not possible in isotropic 
materials (jGoldbergl . llQGOl ). such couplings can manifest themselves for special directions 
in crystals , for in stance, for propagation in the [1 1 1] direction in a cubic crystal (see 
Domanskil f 2000al )). These special directions are the acoustic axes. 



In this paper we clarify what kind of coupling is possible according to the type of 
a symmetry axis and we demonstrate how many constants are needed to describe the 
coupling of pairs of (quasi-)transverse waves in a particular case. It turns out that we need 
two constants for 2-fold symmetry axis, and only one constant for 3-fold axis. Moreover we 
also prove that there cannot be a quadratically nonlinear coupling for shear wave equations 
if the symmetry axis is 4-fold or 6-fold. In the absence of symmetry, the propagation of 
shear waves along an acoustic axis depends upon four constants governing the nonlinear 
terms in the coupled equations. It is possible that a pair of evolution equations decouples 
if the four constants satisfy certain special relations, which are derived here. We illustrate 
these general statements with some examples of particu lar elastic materials. Some of 
these results appeared previously in an abbreviated form (IDomanski and Norrisl . |2008| ). 

The paper is organized as follows. In Section [2] we present the model of nonlinear 
elastodynamics and its constitutive assumptions and we demonstrate the reduction of the 
governing equations to a quasilinear plane wave system. Section [3] contains a presentation 
of the method of Weakly Nonlinear Geometric Optics (WNGO) and its applications to 
quasi-longitudinal and quasi-transverse waves. Special attention is devoted to the degen- 
erate case of pairs of quasi-transverse waves propagating along acoustic axes. Coupled 
evolution equations are derived for two-fold and three-fold symmetry axes. The analysis 
of the coupled evolution equations is discussed in Section H] under conditions of symme- 
try about the propagation direction. Explicit examples which illustrate some of these 
theoretical results are provided for cubic crystals. Section [5] considers the special case of 
quasi-transverse waves propagating along an acoustic axis in the absence of material sym- 
metry, and derives a condition both necessary and sufficient that the evolution equations 
decouple. Some concluding remarks are offered at the end of the paper. 



2 Preliminaries 
2.1 Basic equations 

The equation of motion of a continuum written in material (Lagrangian) coordinates is, 
in the absence of body forces, 

Po^ = DivT. (1) 

Here po is the mass density in the reference configuration, u is the particle displacement, 
T is the first Piola-Kirchhoff stress tensor, and Div denotes the divergence operator with 
respect to the material coordinates X. In the hyperelastic medium there exists a stored 
energy density per unit volume in the reference configuration, denoted by iy(F), such 
that 
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where the deformation gradient is 

F = I + Vu (3) 

with I the identity tensor. 

2.2 First order system 

It is convenient for our purposes to write the equation of motion as a first order system 
of partial differential equations. To this aim we introduce the particle velocity 



du 

and instead of ([H) we write the system of equations 



(4) 



po^=DivT, (5a) 
^ = Vv. (5b) 

Equation ( I5b[) follows from the comparison of the time derivative of (I3l) with the space 
gradient of (jl]), assuming that the displacement vector u is at least twice continuously 
differentiable. We introduce the following notation for the displacement gradient 

M = Vu. (6) 

Then, using ([2]) and ([3]) in the form F = I + M we can express the system ([5]) as 

2.3 Constitutive relations 

Components of vectors and other quantities will be referred to an orthonormal basis 
{ei, 62, 63}. The strain energy in an arbitrary material is assumed to have the following 
expansion for small strains: 

1 1 

^ ~ 2! '^'^'"^'^ ^""^ ^^'^ ~'~ 3! '""^'^'^'^■^ "^"^ ^^'^ + ■ ■ ■ ) (8) 

where Eab are the components of the left Cauchy- Green strain tensor 

E = i(F*F-I) = ^(M + M* + M*M), (9) 

and the summation convention on repeated subscripts a,b, ■ ■ ■ ,g,h is assumed. The 
symmetry of the strain implies the relations 

^abcd ^cdab ^bacdi ^abcdef ^cdabef ^cdefab ^bacdefi (-^0) 
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which imply that there are no more than 21 and 56 independent second and third order 
moduh, respectively. The second order moduli are assumed to be positive definite in the 
sen se that Cg b cdSgbS cd > for all non-zero s = s*. 

iHearmoiil ( 19531 ) provides a complete enumeration of the third-order constants for all 
crystal classes. Equations ([2]), ([3]), ([8]) and IQ together imply that the Piola-Kirchhoff 
stress is 

Tab = Cabcd M,d + ^Nabcdef M^M.f + ^Ngbcdefgh M,dM,fMgh + ■■■ , (11) 

where 

^abcdef ^abcdef ~l~ ^abdf ^ce + Cbfcd ^ae + Cbdef ^ac ■ (12) 

Thurston! fjl984l . eq. (38.5)) gives an expression for the higher order coefficients N abode fgh- 



Note that N abode f 7^ Nbaodef, which implies that the non-symmetry of T is a second-order 
effect. 

For the remainder of the paper we take po = 1 for simplicity. 



2.4 Plane waves 



Plane wave solutions are described by displacement u that depends upon a single com- 
ponent of X, say a; = X ■ n, where n is the direction of propagation. The displacement 
gradient of ([6]) reduces to 

M = m(g)n, (13) 



where 



m 



du 

dx ' 



(14) 



is a displacement gradient vector. Defining the energy function for plane deformation by 



V{rn) = W{rn ® n) 
we can write the system of elastodynamic equations 

d\ d'^V dm. 

dt dm'^ dx ' 
dm d\ 



for plane waves as 



dt dx 

The above system can be expressed in a quasilinear form as 

dw 



dw 
'dt 



A(w] 



dx 



0, 



where 



The 3x3 matrix B is 



w 



Afwl 



B 

1 



B = A + *m 



-Ilmm 



(15) 

(16a) 
(16b) 



(17) 



(19) 
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where A(n), ^'(n) and n(n) are, in components, 

Aqc = CabcdXlbnd, 
^ace = Nabcde/nbUdUf, 
^aceg = NabcdefghnhndnfUh, 

or, in short, 

Bac = CabcdnbUd + N abcdefUbUdU f TUe + ^NabcdefghnbUdUfUh niemg H . (21) 

The positive definite property of the second order moduh imphes that A, also known 
as the Christoffel or acoustical tensor, has spectral form 

A = aiki ® ki + a2k2 ® kg + askg ® kg, (22) 

where aj > and {ki, k.2, kg} is an orthonormal triad of vectors. The six eigenvalues of 
A(0) therefore split into three pairs with opposite signs: 





/«1 = 


-A2, 


(23a) 


A3 = -A 


/tt2 = 


— A4, 


(23b) 


A5 = -A 


/os = 


-Ae. 


(23c) 



The corresponding right and left eigenvectors of A(0) are, respectively, 
^^-^-^ k, J' '^^-[ k, J' 

J = 1,2,3. (24) 

Note that Ij ■ = 5,^. 

Remark 1. The property a^- > implies that all six eigenvalues A^, = 1, 2, . . . , 6 of A(0) 
are real. We assume that all the right and left eigenvectors of A(0) form a full set of 
linearly independent eigenvectors. The system is therefore hyperbolic at the origin. This 
assumption can be expressed in terms of restrictions on the strain energy (e.g. rank-one 
convexity), but we will not discuss this issue here. We would like to emphasize, however, 
that we do admit the possibility of coincident eigenvalues, i.e. non-strict hyperbolicity. 
Moreover, the main objective of this paper is the analysis of the case when pairs of 
coincident eigenvalues correspond to (quasi-)shear waves. We call such waves degenerate. 

3 WNGO asymptotics 

In this section we are interested in deriving the simplest nonlinear evolution equations 
for the amplitudes of weakly nonlinear elastic plane waves. Using the method of weakly 
nonlinear geometric optics (WNGO) we will first derive equations for single waves and 
then for coupled pairs of waves. No attempt is made here to review WNGO, which 



(20a) 
(20b) 
(20c) 
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has been widely deve l oped and applied in many areas of mechanics and physics since t he 
seminal work by iLaxI ( 19571 ). We refer the interested read to the review of iHunted fll995l l 



For single waves two cases will be discussed. In considering the system (1171) we first 
choose an eigenvalue As(w) of the matrix A(w) in (fT7|) such that 

VwA,(w) -r.lw^o ^0, (25) 

where is the eigenvector corresponding to A^. This assumption, called genuine nonlin- 
earity, at the zero constant state is typically satisfied for longitudinal or quasi-longitudinal 
waves. Using the perturbation method we will derive a nonlinear evolution for such waves 
in Sec. 13.11 Next, in Sec. 13.21 we also present the simplified evolution equations for single 
waves which do not satisfy the assumption fl25l) . This is typical for shear or quasi-shear 
elastic waves. Finally, in Sec. 13.31 we will also derive the evolution equations for the 
amplitudes of coupled waves that have coincident speeds. 

3.1 Evolution equations for single (quasi-) longitudinal waves 

Let us consider an initial- value problem for a quasilinear hyperbolic system (fTTIl : 

-W + ^f"- ' ar = »• (26) 

w''(0,x) = e wi(0, X, x/e), 

where e is a small parameter. Note that the high frequency ('fast') variable ^ appears in 
the initial condition. The initial data are presumed to have compact support. 

We consider the single wave weakly nonlinear geometrical optics solution to fl26|) 

w^(t,x) = ea,(t,x,r/)r, + 0(e2), (27) 

with an unknown amplitude and a new independent variable 

= e^^{x - \st). (28) 

It is assumed that the eigenvalue A,; is distinct from the others, and its eigenvector is r^. 

Remark 2. Instead of using a 'fast' variable, we may equally well use a 'slow' variable ex 
and work with a solution of the form w = w{x — At; ex), which represents a traveling wave 
modified by nonlinear effects over a large length scale. Using this alternative starting 
point we could obtain the same results as WNGO. 

Introducing the ansatz (1271) into (!26|) . using a Taylor expansion of A(w^) around 0, 
we then apply the method of multiple-scale asymptotics. This relies on treating 77 as a 
new independent variable. We sequentially collect terms of like powers in e and equating 
these terms to zero one finds that the solvability condition applied at the 0(e) level yi elds 



a nonlinear evolution equation for the unknown amplitude as (see ( iDomanskil . l2000af ) for 
the details): 

where the coefficient of nonlinearity is 

= 1, ■ (VwA(w)r,r,) |^^q. (30) 
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One can show that Tg is equal to the left hand side of (1251) . and hence Tg 7^ 0. An explicit 
equation for the nonlinearity coefficient follows from eqs. (fT2l) and (fT7j) as 

Tg = — A„ "'"kis+ii ■ ^ k|s+i|k|s+i| 

1 I ^zLij I s+lj I s-|-lj 2 

= -X'^CahcdefUbUdUfka ^ kc ^ ke ^ + -A^n ■ k[^s±lj , (31) 

where [ J denotes the floor function, that is, the largest integer less than or equal to a 
given number. We also use the notation k^ = /cjei + ^2^2 + ^3^3 in the Cartesian basis. 
The parameter Tg is well known in n onlinear elast odynamics - for instance, it governs the 
growth of elastic acceleration waves flNorrisl . [l99lh . 

Example 1 . Consider an isotropic elastic medium for which the strain energy is given by 
W = + 2fiL)ll - 2fiLlI^ + ^{Im + 2mM)ll - ^niMUh + umIII^, (32) 

where XliJ^l are second order Lame constants, Im, '^m and um are third order Murnaghan 
constant^], and the strain invariants are 

/e = tr E, //e = ]- [(tr E)2 - tr E^] , III^ = det E. (33) 

In this case the two coefficients appearing in the evolution equation fl29|l are 

, / , , ^ 3(Al + 2fiL) + 2{Im + 2mM) 

A, = - V Al + 2/iL and Tg = ■. . (34) 

Hence, we see that in the isotropic elastic medium the coefficients in the evolution equation 
for longitudinal waves are determined by the two second order and the two third order 
constants, Xl, /ii and Im, itT'M, respectively. 



3.2 Evolution equations for single (quasi-)shear waves 

We now present the simplest nonlinear evolution equation for the amplitude ag of a single 
weakly nonlinear (quasi-)shear elastic wave for which Tg = 0. The case when = 
requires a different scaling from that of (l28l) , one appro priate to cubic non linearity as the 
leading term. The procedure is described in detail by iDomanskil f 2000al ). The WNGO 
solution for the single wave has the same formal expansion as in eq. (1271) . but now, 
crucially, 77 = e~^{x — Xgt). The modified asymptotics leads to the following governing 
equation for ag-. 

d(jg , dcjg 1 _ dal „ 

(35) 



, dcTg 
+ Xg 



-G ^ 
3 drj 



0, 



dt dx 
with nonlinearity coefficient 

Gs = ^X;\s±,^ ■ (3*kL^jq + nkL^jkLs±ijkLs±ij), 

where the vector q is orthogonal to k^ and satisfies 

(A-A2l)q + *kL^jkL^j =0. 



(36) 



(37) 



^Relations between the t hree Murnaghan con stants and alternative triads of third order constants for 
isotropic soUds are hsted bv lKostek et al. (1993). 
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Example 2. For an isotropic elastic medium with strain energy given by (!32|) the two 
coefficients in the evolution equation (l35l) are 



A. 



/iL and Gs = — 
4As 



>^L fJ'L + {fJ'L + 



(38) 



Hence, in the isotropic case, the coefficients in the evolution equation for shear waves are 
determined by the two second order and only one of the third order constants, A^, fii, 
and rriM, respectively. 

Remark 3. By applying the method of characteristics we reduce the differentiation -ir+XS- 



to the differentiation where r 

OT 

partial differential equation 



dt dx 

Xx is a characteristic variable. In this way the 



dt dx j drj 



(39) 



for the function of three independent variables a 
differential equation 

da 



l^da^ 



j drj 



a{t,x,ri) transforms to the partial 



0, (40) 



for the function of two independent variables a = a(r, r]). Here j = 1, 2, 3, ... is a natural 
number. When j = 2, as occurs for the longitudinal wave evolution equation 0291) . we 
call the equation ( HOll the inviscid Burgers equation. Similarly, j = 3 for the (quasi-) shear 
wave in (1351) . the equation (HUj) is called the modified inviscid Burgers equation. 



3.3 Degenerate plane waves: acoustic axes 

In this section we derive the simplest nonlinear evolution equations for the amplitudes 
of a pair of weakly nonlinear elastic waves in the case when these waves have coincident 
wave speeds. 

Consider the Christoffel tensor A from fl22l) . 
Definition: We say that the eigenvalues of A are degenerate, if 

A = ai (I - kg ® kg) + agkg ® kg, (41) 

that is, if a-\ = a '?, i n fl22l). In such a situation we say that n is an acoustic axis (see 
KhatkevichI f!l962l ) or lNOTrisI (120041 )). 



Let {ki, k2, kg} be an orthonormal triad of vectors and define the left and right 
eigenvectors of A(0) as before, see eq. (l24l) . Then A(0) has two pairs of coincident 
eigenvalues: A^ = Xs+2 for s = 1 and s = 2. We now consider the following ansatz for the 
initial value problem (!26|) : 

w{t,x) = e{asit,x,r])rs + as+2it,x,r])rs+2) +0{e'^), s = 1,2, (42) 



with 7] = e ^{x — Xst), and where as and as+2 are the unknown amplitudes. Inserting 
(142|) into (!26|) and applying the multiple scale asymptotic methods described in Section 
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lO (see [D omanskil (j2000al )) we obtain the following pair of coupled evolution equations 
for the amplitudes ag and (Ts+2- 

+ A.^ + V r^,.^ + 2riJ-^ + ru.^sJ4^) - o, (43a) 



dt ' dx 2\ dr] dr] dr] 

^ + + ^ f r:f ^ + 2rn^.^^^ + r-^ .J4^) = o, (43b) 

dt dx 2\ dr] *''+^ dr] s+2,s+2 j \ j 

where the interaction coefficients are, in general, 

r^,, = 1. ■ (Vw A(w)rpr,) |^^o- (44) 

The assumption of hyperelasticity ([2]) implies that these coefficients have the following 
symmetry property 

n,, = ^ip- (45) 

Moreover the formulas (12^ imply that = ~^Pq^ J = 1,3,5. In our case we have 

Using Cartesian components we can express the interaction coefficients as follows 
1 

2' 



1 _-|^ p [i+ij [£±ij [2±ij 



[^i+ij [^E+ij [^£±ij 

"I" CabcdXI'bf^cf^dika ^ ^pg -\~ ka ^ Sjg + /Ca ^ ^jpjj 



1 _1 1 — 1 |2±i| |2±i| 

= -Aj CabcdefnbridUfka ' kc^ fee ' 

+ ^Aj-n ■ (k^i+ij(5pg + k^p+ijjjg + k^a^^Sjp) . (47) 

Therefore the coefficients have in our case, in addition to the general property the 
following indicial symmetries: 

This makes them totally symmetric under the interchange of indices. In particular, note 
that 

r^+' = = ^K-'cabcdefnbndnfk![^klr^kl[^^ + ^A,n ■ k^^j = (49a) 

r:+2,,+2 = r:+^2 = IXs^'CabcdefribUdUfk}^^ kl^^^ kP^^ + ' ^^L^J = ^^+2" (49^) 

In summary: 

Lemma 1. The evolution equations for the amplitudes of shear waves propagating along 
an acoustic axis are 

r:«^) = 0, (50a) 

The nonlinear terms in the equations involve four coefficients: Tg and Ts^2 from eq. I[3l\) . 
r^"*"^ and T'l^2 from eqs. (^^p, the latter two of which determine the coupling between the 
equations. 



dx 


4l 




, o-ps+2^(^«^^+2) 
dT] 


das+2 1 




.,+2d(r'i 


, o-p« d{(rsCrs+2) 
^ dr] 


dx 




' dr] 
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4 Simplification along symmetry axes 

The previous analysis shows that there are at most four independent coefficients appearing 
in the nonhnear terms of the shear elastic wave's coupled system fl50l) . We will demonstrate 
the following: 

Lemma 2. The number r of coefficients in the nonlinear terms in the coupled equations 
( f5^) for a pair of shear waves is always reduced if the direction of propagation n is a 
symmetry acoustic axis. Specifically, 

r = 2 for propagation along a 2-fold symmetry acoustic axis: Tg = rf_,_2 = 0; 
r = 1 for a 3-fold symmetry acoustic axis: Ts = r*_,_2 = and r^,+2 + r^"''^ = 0; 
r = for an acoustic axis of 4-fold or higher symmetry: Tg = r^_,_2 = Ts+2 = 

2 = 0. 



4.1 Two-fold axis 

We say that the propagation direction is a two-fold axis of symmetry if it lies in a plane of 
symmetry of a monoclinic solid. Let e be the normal to the plane of monoclinic symmetry, 
then e also belongs to the plane of degenerate wave vectors. Let e coincide with ei, and 



be parallel to k|^s+ij . It then follows that n ■ k|^s+ij = 0, and from eqs. fl3T|) and (I49b|l . we 
have 

= r^^, = ^Xs-'cabcdefnbndnfk]^^k['^^kfl^\ (51a) 
s s -'- -1 L^^J L""*"^] [""'"^j 

— '^s+2,s+2 ~ 2^'^ Cabcdef^bf^df^fkckeka ^ kc ^ ke ^ ■ (51b) 

Each term in these expressions involves an element of Cabcdef with the index 1 occurring 
either once or thrice. But by definition of a plane of symmetry these elements vanish, and 
hence the two elements and F*_,_2 vanish. The canonical form of the evolution equations 
is 

4.2 Three-fold axis 

The generic configuration in this case is propagation along the axis of a trigonal material. 
This possesses three planes of symmetry each containing the axis of symmetry, and mu- 
tually at 120° to one another. Reasoning similarly to the case of the two-fold symmetry 
axis, we get that F^ = F^^2 = 0; but here we have moreover that 

F,+2 + r:+2 ^ 0, (53) 

see below for details. Therefore the canonical form of the evolution equations for a three- 
fold symmetry axis is 

das , . das -p 9(cT,a,+2) ^ . 

-KT + r,+2 7^ = 0, (54a) 

ot ox orj 
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Hence only one coefficient Ts+2 characterizes the nonhnear terms in the evolution 
equations for a pair of degenerate plane waves in the case when the propagation direction 
is a three-fold symmetry axis. 

Remark 4. The system ( IMI) can be transformed into th e single comp lex Burgers equation. 
This equation was studied in the Noelle's Ph.D. thesis flNoelld . [1990 ) . It was shown there 
that the solutions of the complex Burgers equation become singular in finite time for a 
large class of initial data. The shocks which develop are of nonclassical type, that is, they 
do not satisfy Lax's conditions. This is due to the fact that too few characteristics enter 
the shock front in comparison to compressive shocks which satisfy Lax's condition. 



4.2.1 Proof of equation (dl) 

The basis for the proof rests on the property fHHjl combined with the three-fold symmetry. 
Using the same arguments as for the two-fold axis, the three-fold symmetry is associated 
with three planes of monoclinic symmetry with normals all perpendicular to n and since 
each normal lies in the plane of degenerate wave vectors, it follows that the degenerate 
wave vectors are orthogonal to the axis and are therefore pure transverse waves. 

In order to simplify matters, recall that for propagation along an acoustic axis s = 1 
or s = 2 (see eq. (H2l) ). and therefore the coefficients F^^ may be identified with the 
elements of a totally symmetric third order tensor of dimension 2. Thus, 

g = gafi-yK ® k^ (g) k^, (55) 

where lower case Greek subscripts take values 1 and 2, and, for instance, 

gap-y = CabcdefnbndUfk'^k^kJ + Cahcdnbncnd{k'^5p^ + k^d^^ + k^Safi) (56) 

In this case the relationship with the nonlinearity coefficients ^ is defined by 

2l;^L4^JL^^JL^J' ^^^^ 

where j,p,q take values 1,2,3,4. The totally symmetric property means that Qaf^^y is un- 
changed under any permutation of the three indices, and hence there are at most four 
independent elements. Consider the change of basis 

k.\ = cos^ki + sin^k2, kg = — sin^ki + cos^k2, (58) 

and define 

9apy{G) = gk'ak^k^, (59) 

then 





3c^s 3cs^ 




(9iiM\ 


9222{0) 


—s^ 3cs^ — 3c^s 




9222(0) 


giuie) 


—c^s cs^ c" — 2cs^ 2c^s — 




^7112(0) 


\9i22{e)/ 


cs^ (?s — 2(?s (? — 2cs'^ 




\9122{0)J 



(60) 



where c = cos 6', s = sin 9. We note the property QaP'^iO) = — ^'0/37(6' + vr) for each element 

Let ei be the normal to one of the three planes, then the coefficients 5'iii(0) and 
5'i22(0) vanish. Similarly, (^iii(±|7r) and (7122 (±|vr) must vanish. Using ( 160|) with5fiii(0) = 
5'i22(0) = 0, we find that 



giui±l7v) = 3(7i22(±i7r) = ±^(^7112(0) +^?222(0)), 



(61) 
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and hence the 3-fold symmetry requires that 

^7ll2(0)+^?222(0) = 0. (62) 

This is precisely eq. f l53|) . In summary, the evolution equations remain coupled, but are 
characterized by a single nonlinearity parameter and take the form fl5^ . 



4.3 Four- fold axis 

If the axis is one of 4-fold symmetry, then the same arguments as above imply that the 
degenerate wave-vectors are orthogonal to the axis, and are thus pure transverse waves. 
The same arguments also imply that the elements of F associated with the second plane 
of symmetry must also vanish, and hence all four are zero. Thus, F = and there is no 
coupling. 

The same reasoning applies to axes of higher symmetry, since all such axes are equiv- 
alent to an axis of transverse isotropy. Therefore the canonical form of the evolution 
equations for a symmetry axis with four-fold or higher symmetry is 

^ + = 0, (63a) 

dcrs+2 . . das+2 ^ /«qu\ 

4.4 Applications to cubic crystals 

Let us consider a cubic crystal of class m3m in which th e strain en e rgy W is defined by 
three second order and six third order elastic constants ( Domanskil . l2000al lbl. l2008l ). (see 
also the Appendix): 

W = W{cn, ci2, C44, Cm, C112, C144, C123, ciee, £455). (64) 

We focus on three particular directions of plane waves' propagation: [10 0], [110] and 
[111]. 

Example 3. Consider first the case of propagation direction n = [10 0]. This direction is 
along a four-fold symmetry acoustic axis. The shear wave speeds are given by 

Ai = A3 = —y/c^ = — A2 = — A4. (65) 

Hence, the eigenvalues Ai, A2, A3 and A4 are coupled but the second order nonlinearity 
coefficient Tg vanishes for s = 1,2,3,4. Therefore there is no quadratically nonlinear 
coupling in the evolution equations for shear waves. One can show that propagation of 
the s hear plane waves is described by the evolution equation with cubic nonlinearity, eq. 
f l35|) fjPomanskil . l2000a ). The coefficients at cubic nonlinearity in this equation are as 



follows: 

2(cii + ciee)^ + 044(011 + ciee) + Ci66(c44 + ciee) 



4A. 



Cii — C44 



5 = 1,2,3,4. (66) 



Example 4. Consider now the case of propagation direction n = -^[110]. This direction 

is along a two-fold symmetry axis and is not an acoustic axis. The shear wave speeds are 
given by 

Ai = -\r^^^-Tr^ = -A2, A3 = -y/c^ = -A4. (67) 



Degenerate nonlinear elastic waves 



December 18, 2008 



13 



Hence, the eigenvalues Ai, A2, A3 and A4 are distinct and the second order nonhnearity 
coefficient Tg vanishes for s = 1,2,3,4. There is no quadratic couphng between the 
shear waves, and furthermore one can show that propagation of the s hear plane waves i s 
described by the evolution equation with cubic nonlinearity, eq. (135!) (jPomanskil . l2000a ). 
The coefficients in this equation are: 



16 A. 



[Cll - C12 + |(ciii - Cii2)]^ + 2(cii - Ci2)(ci2 + C44) 



C12 + C44 



s = l,2- (68) 



16A.. 



(Cii — C12 + C144 + C166 + 2C456)^ + 2(Cii — Ci2)(Ci2 + C44) 



3,4. 



C12 + C44 

(69) 

Example 5. Finally let us consider the case of propagation along an axis of three-fold 
symmetry, n = [1 1 1] , which is also an acoustic axis. In this case the shear wave 

amplitudes are coupled and are described by the coupled evolution equations with 
(see the Appendix) 



Ai = A3 = — — ^ ^^'^ = -A2 = -A4, (70a) 

r, = Ti = = 12, (70b) 

= Ta = /- [ciii + 2ci23 ~ 2c456 ~ 3(cii2 — Ci44 + Cige)] = — C^'Oc) 
lov 2Ai 

5 Degenerate transverse waves in the absence of sym- 
metry 

5.1 The principal result 

In this Section we consider the general case of quasi-transverse degenerate wave vectors. 
No symmetry is assumed. We examine the possibility that the coupled nonlinear wave 
equations for the two amplitudes decouple, and derive a general condition that is both 
necessary and sufficient for this to occur. 
The main result is the following: 

Lemma 3. A pair of quasi-transverse waves propagating along an acoustic axis are de- 
coupled if and only if the nonlinearity coefficients satisfy the identity 

T,TU2 + r.+2rr^ - (r:+2)' - (r^')' = o. (7i) 

// this condition is met then there is a coordinate transformation for which the coupling 
terms disappear and the degenerate transverse waves satisfy separate but different evolution 
equations: 

da. , da. 1 ^, da'^ ,^ , 



da^^o , 9cr.,,9 1 , da 



2 



ir + A.^+2^:«-^-0, (72b) 
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Remark 5. The condition (I7T1) can be expressed in terms of the 2- dimensional third order 
symmetric tensor g as (see eq. (l84l) ) 

galSygaP^ — galSpgaff = 0. (73) 

Based upon eq. fl56p this may be interpreted as a specific constraint on the third order 
moduh involving the direction n of the acoustic axis. It also depends upon the second 
order moduli through the common eigenvalue, but this dependence disappears when the 
degenerate waves are purely transverse, in which case (ITTi) is strictly a relation between 
the third order moduli. 



5.2 Derivation of Lemma [3] 

We use the property seen previously in Section (jlj) that the coefficients T^^ ^ are the ele- 
ments of a third order totally symmetric tensor in two dimensions, g of eq. (l55ll . Let 
Ai = A2 = A. The coupled wave equations can be expressed in terms of a 2- vector for the 
displacement gradient vector of eq. (fT4l) : m = (Tsk^ + as^2^s+2, or m = miCi + 771262, 
where ei, 62 are ortho normal in the plane of k^, \<.s+2'- 

^ ^^~dx' ^ 2 ^7/ ^^^""^^ ^ 25(11277717712 + 5l22"^2) = 0' 

dm2 dm2 1 d , 2 2\ 

+ -^2-^ + 2 ^ '^5'112"^1 + 25122 77117772 + 5222"^2) = 0) 

or 

dm , dm Id, \ ^ 

The coupling between 77^1 and 77^2 vanishes iff the two elements 5112 and 5122 are simulta- 
neously zero. Even if these are non-zero there might exist a coordinate transformation in 
which the transformed quantities vanish. It is this possibility that we seek. 

Thus, we consider the possibility that there is some angle of rotation 9 such that 
gii2{0) = and 5122(6*) = where gap-^id) are defined by (l59ll . These two conditions are 
simultaneously satisfied if the 2- vector ge^eg vanishes. Under the rotation ( |58l) we have 

ge'ie2 = a cos 2^ + b sin 29, (75) 

where 

a = geie2, b = -(ge2e2 - geiei). (76) 

The form of (1751) indicates that the vector ge'^e'g can be zero if and only if a and b are 
parallel. Thus, the general condition for no coupling is 

no coupling ^ a x b = 0, (77) 

which is equivalent to 

no coupling ^ /i = 5^12 + 5122 - 5'ii2fi'222 - 5'i225'iii = 0. (78) 



This is precisely the result of eq. (17T!) in Lemma [31 We will show next that the quantity 
yU of eq. (!78|) is an invariant, independent of the coordinates used. 
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5.3 Tensor properties of the nonlinearity coefficients 

The third order tensor g is 2-dimensional and totally symme tric, that is, the C artesian 
components are unchanged under permutation of the indices. iJerphagnon ( 1970l ) consid- 
ered the general form of third order tensors in 3-dimensions, and based on these results 
we may partition g as follow^ 

g = g{i)+g(3)^ (79) 

where ^ 

g*^^-* is pseudovector (jjerphagnonl . Il970l ) with two independent components, ((7122 + diu) 
and {gu2 + 9222)- g^^\ which also has two independent elements, (8(7122 — fi'iii) and 
(3(7112 — 9222), may be called a dimer. More importantly, it is harmonic, i.e. g^^^^ = 0. 
Hence, the quadratic invariant gap-yQaP') can be expressed 

ga^^ga^^ = 71+73, (81) 

where 71 and 72 are the quadratic invariants of the constituent tensors: 



71 = 9aP^9aP'y = ^(^111 + ^122)^ + ^(^112 + 9222^, (82) 

73 = gSi^gSi^ = ^ (3^122 - ^m)^ + ^ (3^112 - ^222) ^ (83) 



Comparing these with eq. fl75]) indicates that fi = ^73 — |7i, and hence the no coupling 
condition can be expressed in terms of invariants as 71 = 373. Alternatively, noting that 

7i = |fl'Q/3/3fi'a77, we deducc /i = ^ap^gap^ " ^ga^/sga^^, and hence 

no coupling ^ gap^gap-i = gaf3l3gayy (84) 



Example 6. We check the condition eq. fl7T]) in two cases of acoustic axes: [10 0] and 
[1 1 1] for the cubic crystal considered earlier. It is easy to see that the condition fl7T]) is 
satisfied and the evolution equations for shear waves are decoupled for the [1 0] acoustic 
axis. However, for the [1 1 1] axis, the left hand side of eq. flTTj) is equal to — 2rs_|_2 ^ 0, 
so the condition fl7T]) is not satisfied. Therefore there is quadratically nonlinear coupling 
in the evolution equations for shear waves in this case (see eqs. fl5^ and fITOl) ). 



6 Concluding remarks 

Starting from a formulation of the governing equations as a first order system of quasi- 
linear equations, we have derived the general form of the amplitude evolution equations 
for weakly nonlinear plane wave propagation. The major new results concern the form of 
the evolution equations for the degenerate conditions associated with propagation along 
acoustic axes, summarized in Lemma 1 and eqs. (l50l) . The quasi-transverse wave am- 
plitudes are coupled at the quadratically nonlinear level, with at most four interaction 
coefficients. The number of coefficients reduces in the presence of symmetry, with the 



^Ge neral third carder t ensors include additional terms g*-"-* and g*-^-* which vanish for a symmetric 
tensor ( Jerphagnon . 1970l ). 
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precise number determined by Lemma 2. For instance, the coupling in the presence of 
three-fold symmetry about the propagation direction depends on a single interaction coef- 
ficient, with the canonical form of the coupled equations given by eqs. (154|) . The nonlinear 
coupling disappears if the acoustic axis has four-fold or higher symmetry. The isotropic 
solid is of course the most obvious example, but the results presented here show that 
similar decoupling can be expected in the presence of anisotropy. We have also shown 
that it is possible for the coupling to vanish even when the acoustic axis is not a symmetry 
axis. The condition, defined by Lemma 3, requires that the interaction coefficients satisfy 
a unique relation. Taken together, the variety of results presented here shed light on the 
nature of the equations governing nonlinear wave propagation in elastic crystals. 
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A Appendix: Cubic crystals 



In a cubic crystal with cube axes e^, j = 1, 2, 3, we have (see iDomanskil tOO^ ) 

CabcdEabEcd = Cll(/i - 2I2) + 2C12/2 + 4C44/3, 
CabcdefEabEcdEef = Ciii(/f — 3/i/2 + 3/4) + 3Cii2(/l/2 — 3/4) 

+ 12Ci44(/i/3 - h) + 6C123/4 + 48C456/5 + 12Ci66/6, 



where E = E and 

h 
h 
h 
h 
h 
h 



Ell + E22 + -£"33, 

E11E22 + E22ES3 + i?33-E'ii, 

E12 + ^23 + E^i, 

E11E22E33, 

E12E23E3I, 

{Ell + E22)El2 + (-E'22 + E33)E23 + (£"33 + Eii)El^. 



The four cube diagonals are axes of trigonal or 3-fold symmetry, and are acoustic 
axes. Consider propagation along the cube diagonal acoustic axis n = [1,1,1]/V3, and 
assume s = 1. The parameters Tg and rs_|_2 follow by taking ki orthogonal to one of the 
3 symmetry planes, e.g. ki = [1,-1, 0]/v^ and k2 = [— 1, — 1, 2]/-\/6. The coefficient Fi 
is then obtained by using 



E= -(n®ki + ki®n) = ^ 



2 

-2 

1 -1 



1 



Therefore, h = h = I5 = Iq = implying CabcdefEabEcdEef = and hence F^ = 0, as 
expected for a two-fold symmetry axis. The coefficient F3 = Ts+2 follows from eq. (13T1) 
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with 



E = ^(n ® k2 + ka O n) = 

2 Qy/2 



-2 -2 1 
-2 -2 1 
1 1 4 



leading to 



CabcdefEabEcdEef — [cm + 2Ci23 — '^C^^Q — 3(cii2 — C144 + Ciee)] • 
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